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ABSTRACT 

We analyze stellar tidal disruption events as a possible observational signature of 
gravitational wave induced recoil of supermassive black holes. As a black hole wanders 
through its galaxy, it will tidally disrupt bound and unbound stars at rates potentially 
observable by upcoming optical transient surveys. To quantify these rates, we explore 
a broad range of host galaxy and black hole kick parameters. We find that emission 
from a transient accretion disk can produce ~ 1 event per year which LSST would 
identify as spatially offset, while super-Eddington tidal flares, if they exist, are likely 
to produce ~ 10 spatially offset events per year. A majority of tidal disruption flares, 
and a large majority of flares with an observable spatial offset, are due to bound rather 
than unbound stars. The total number of disruption events due to recoiled black holes 
could be almost 1% of the total stellar tidal disruption rate. 
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1 INTRODUCTION 

Recent advances in numerical general relativity quantified 
how the coalescence of unequal black hole binaries leads to 
the anisotropic emission of gravitational radiation, which 
can carry enough linear momentum to deliver a substantial 
kick to the merged black hole (Pretorius 2005, Baker et al. 
2006, Campanelli et al. 2006). Because the inspiral and even- 
' tual coalescence of supermassive black holes (SMBHs) is a 
$H ■ frequent consequence of galaxy mergers, it is expected that 
SMBH recoil will often accompany these events. The magni- 
tude of the recoil varies as a function of the initial mass ratio 
and the spin and orbital angular momentum vectors of the 
coalescing black holes, but is generally of order hundreds 
of kilometers per second (Schnittman & Buonanno 2008, 
Lousto et al. 2010). At the high end of the velocity distribu- 
tion, the black hole can recoil with a velocity > 1000 km s~^ 
and escape the merged galaxy. More commonly, it will os- 
cillate for a time ranging from 10® - 10^ yrs (Blecha & Loeb 
2008, Sijacki et al. 2010, Blecha et al. 2010) before settling 
down into the galaxy's center. 

Observation of a recoiling SMBH would be of interest 
both as a probe of hierarchical galaxy evolution and as a test 
of the predictions of strong field general relativity. Conse- 
quently, several papers have focused on observational signa- 
tures of a recoiling SMBH. A black hole kicked at a substan- 
tial fraction of its host galaxy's escape velocity can create 
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or expand a low-density stellar core (Gualandris & Merritt, 
2008). The small cloud of stars gravitationally bound to the 
recoiling black hole would appear like a globular cluster, but 
with much higher velocity dispersion (O'Leary & Loeb 2009, 
Merritt et al. 2009, O'Leary & Loeb 2011). Gas accretion 
onto the black hole, manifested as a spatially or kinemat- 
ically offset quasar (Madau & Quataert 2004, Loeb 2007), 
is potentially a very clear signature, but the gas reservoir 
bound to the black hole will be depleted within 10^ years 
after ejection from the center of its galaxy (Blecha & Loeb 
2008, Blecha et al. 2010). Several candidates have already 
been detected with this last method (Komossa et al. 2008, 
Shields et al. 2009, Civano et al. 2010), although the ev- 
idence for black hole recoil is not conclusive (Bogdanovic 
et al. 2009). A potentially longer-lived source of accretion 
power is the tidally disrupted matter from stars passing too 
close to the recoiling black hole, which could be visible as 
an off-center tidal disruption fiare. Because tidal disruption 
flare lightcurves have, in principle, several unique identifiers 
(Strubbe & Quataert 2009, Guillochon et al. 2009, Kasen 
& Ramirez- Ruiz 2010, Strubbe & Quataert 2010), an off- 
center or intergalactic tidal disruption fiare would be a very 
strong indication of a recoiling black hole. This scenario was 
first investigated in a paper by Komossa & Merritt (2008, 
hereafter KM08). However, to evaluate the utility of tidal 
disruption signatures for recoiling black holes it is first nec- 
essary to obtain estimates for the tidal disruption rates pro- 
duced by these systems. KM08 calculate these rates for rela- 
tively large kick velocities, while in this paper we generalize 
the calculation to cover realistic, and often moderate, kick 
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velocity distributions. The aim of these calculations is to ap- 
proximate the frequency of off-center and spectrally shifted 
tidal disruption flares, to determine if tidal disruption events 
(TDEs) are useful probes of physically motivated kick veloc- 
ity distributions. These rate calculations are also relevant to 
SMBHs ejected from galaxy centers in triple-SMBH interac- 
tions (Hoffman & Loeb 2007), although we do not attempt 
to evaluate the more complicated distributions of parame- 
ters for such events. 

Very large increases in tidal disruption rates have been 
discussed previously in the context of galaxy mergers, both 
as a prompt electromagnetic signal immediately following 
SMBH coalescence (Stone & Loeb 2010), and due to reso- 
nances or chaotic stellar orbits at the end of the dynamical 
friction phase of an SMBH binary (Ivanov et al. 2005, Chen 
et al. 2009, Wegg & Bode 2010, Chen et al. 2011). Neither of 
these mechanisms is capable of producing a tidal disruption 
flare with an observable spatial offset, as both can only oc- 
cur in galactic nuclei. Tidal flares produced by wider SMBH 
binaries could still be a source of confusion in interpret- 
ing spatially offset TDEs, however. Binary-produced flares 
would themselves be of interest, but would not be directly 
useful for testing predictions of gravitational wave recoil. 
Distinguishing between the two possibilities will be simple 
in some cases: for instance, a TDE in intergalactic space, in- 
dependent of an observable stellar population, would likely 
be due to a recoiling black hole; whereas a TDE with peri- 
odic interruptions in its lightcurve (Liu et al. 2009) would 
be due to a hard SMBH binary. However, in many cases 
it will be not be trivial to disentangle the two causes, and 
detailed observation and modeling of the host galaxy could 
be necessary to determine if it is likely to harbor a binary 
SMBH at the observed spatial or kinematic offset. A final 
source of confusion could be disruption of stars by IMBHs 
left over in a galaxy's halo from earlier stages of hierarchi- 
cal growth (O'Leary & Loeb 2011). However, these TDEs 
would probably be distinguishable due to their low black 
hole masses, and in any event would be a complementary 
example of black hole recoil. 

Although the current sample of observed TDEs is small, 
with roughly a dozen strong candidates (Gezari et al. 2009), 
current and upcoming time-domain optical transient sur- 
veys such as Pan-STARR^I, PTFQ, and LSSlB are ex- 
pected to increase that sample by 1 — 3 orders of magnitude 
(Strubbe & Quataert 2009), making it profitable to study 
sub-populations of TDEs. If sufficiently large, the subset of 
disruption flares associated with recoiling black holes could 
be used to constrain the LISA0 event rate or the distribution 
of kick velocities associated with astrophysical SMBH merg- 
ers; at the very least, detection of this subset of TDEs could 
qualitatively confirm recent numerical relativity predictions. 

The outline of this paper is as follows. In §2 we de- 
velop the model used to estimate the TDE rate over a wide 
range of kick velocities and galaxy parameters, and in §3 we 
explain in more detail the distributions of disruption prop- 
erties we integrate over. In §4 we discuss the results of our 
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modeling, and in §5 we offer our conclusions on the viabil- 
ity of tidal disruption flares as a technique for identifying 
recoiled SMBHs. 



2 MODEL 

The three primary quantities to be calculated are the tra- 
jectory of the kicked black hole through its host galaxy, the 
rate at which it disrupts unbound stars that it encounters, 
and the rate at which it depletes its cloud of gravitationally 
bound stars via tidal disruption. The location (and possi- 
ble kinematic offset) of a disruption event relative to its 
host galaxy will determine whether the flare can be distin- 
guished from a TDE due to a stationary black hole. §2.1 
examines expected ranges of black hole kicks, §2.2 discusses 
the relevant tidal disruption physics, §2.3 lays out the galaxy 
parametrization used in this paper, and §2.4 and §2.5 de- 
scribe the techniques used to estimate TDE rates for un- 
bound and bound stars, respectively. Finally, in §2.6 we dis- 
cuss observational constraints. Our general strategy is as fol- 
lows: Strubbe & Quataert 2009 (hereafter SQ09) calculate 
the number of TDEs a generic survey would be expected to 
detect, binned by black hole mass (see SQ09 Figures 9, 13). 
In this calculation they assume a time-averaged tidal disrup- 
tion rate N = 10~^ yr~^ galaxy"^. Our model averages over 
a kick velocity distribution and a range of galaxy parame- 
ters (detailed in §2.3) to calculate a dependent on black 
hole mass, and then modifles the SQ09 survey calculations 
accordingly. 

2.1 Black Hole Kicks 

The gravitational recoil velocity of a post-merger SMBH de- 
pends only on the mass ratio, spin amplitudes, and spin ori- 
entations (relative to the orbital angular momentum plane) 
of the two progenitor SMBHs. The merger of Schwarzschild 
black holes represents a simple case, with a maximum kick 
velocity of ~175 km s~^ occurring at a mass ratio of 0.36 
(Gonzalez et al. 2007). As the dimensionless spin parameters 
of the two SMBHs ai, a2 increase, so does the maximum kick 
velocity iik, although it has complicated functional depen- 
dences on the relative inclination of the pre-merger spin axes 
and the orbital angular momentum axis. An exact calcula- 
tion requires the full framework of numerical relativity, but 
analytic fitting formulae can be calibrated to numerical rel- 
ativity results, recently yielding quite high accuracy (Lousto 
et al. 2010). Although certain combinations of SMBH binary 
initial parameters can produce kicks up to ~ 4000 km s~^ 
(Campanelli et al. 2006), observation constrains how fre- 
quent such high-velocity kicks can be, due to the fact that 
galaxies with bulges all seem to possess central SMBHs (Fer- 
rarese & Ford 2005, Blecha et al. 2010). This constraint is 
weakened, however, by the hierarchical nature of structure 
formation (Schnittman 2007). 

Higher velocity kicks become more important if one as- 
sumes high values of pre-merger spin amplitudes: Lousto et 
al. (2009) found 23% of their Monte Carlo sample to ex- 
ceed 1000 km s~^ by assuming all pre-merger spin ampli- 
tudes ai = a2 = 0.97 (with isotropic distribution of spin 
angles), and mass ratios q between 1 and 1/16. Spins of 
this magnitude are likely realistic for a nontrivial fraction 
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of astrophysical black holes: the Fe Ka line has implied a 
near-max;imally spinning SMBH candidate (Brenneman & 
Reynolds 2006). However, there is reason to believe that the 
spins of merging black holes may align during the inspiral 
phase (Bogdanovic ct al. 2007) in the case of gas-rich merg- 
ers. This has the potential to strongly suppress high-velocity 
kicks relative to an isotropic spin distribution (but see King 
et al. 2005 and Lodato & Pringle 2006 for a description 
of how counter-alignment, which does not suppress as 
strongly, can also occur). Dotti ct al. (2009) find a dramatic 
reduction, with median kick velocities below 70 km s~^, al- 
though their scenario investigates the specific case where q 
is very near 1, while wc arc concerned with a wider range of 
mass ratios. Their results also indicate that cooler gas accre- 
tion (adiabatic index 7 = 5/3) more effectively aligns spins 
and suppresses kicks than a warmer (7 = 7/5) accretion 
flow. More recently, it was discovered that even in gas-free 
mergers, relativistic spin precession is capable of aligning or 
anti-aligning progenitor SMBH spin vectors, and reducing 
median kick velocities (Kesden et al. 2010). Consideration 
of all these factors highlights the need to generalize the work 
of KM08 to a more realistic distribution of kick velocities, as 
recoil velocities in excess of galactic escape speed are likely 
to be uncommon even if SMBH spins arc near-maximal and 
unaligned. 

Because the real distribution of kick velocities depends 
on the SMBH spin distribution and the uncertain physics of 
the last stages of SMBH merger, wo calculate black hole tra- 
jectories given kicks of 100, 200, 300, 400, 500, 600, 700, 800, 
and 900 km s~^, and then interpolate quantities of interest 
such as wandering lifetime and time-averaged tidal disrup- 
tion rates. We then fit these functions of Uk to plausible black 
hole kick distributions. In particular, wo consider the three 
distributions in Figure 2 of Lousto et al. 2010, which as- 
sume isotropic pre-merger spins, spins aligned to within 30 
degrees, and spins aligned to within 10 degrees (dry merg- 
ers, hot wet mergers, and cold wet mergers, respectively). 
These scenarios assume a q uniformly sampled between 
and 1. The dry merger scenario assumes spin magnitudes 
randomly sampled between and 0.9, while the wet mergers 
assume spin magnitudes between 0.3 and 0.9. In reality, the 
distribution of spin amplitudes also depends on factors not 
considered; chiefly, relativistic spin precession (Kesden et al. 
2010) and chaotic versus standard accretion (Berti & Volon- 
teri 2008). Likewise, modeling SMBH merger histories and 
spin evolution using the Press-Sclicchtcr formalism indicates 
a top-heavy spin distribution and a bottom heavy q distri- 
bution (Volonteri et al 2004). We leave more complicated 
oi, 02, and q distributions to future research, although the 
kick distributions we consider should bracket a wide range 
of the available parameter space. 

2.2 Tidal Disruption Physics 



Stars that pass within a radius 



n 



(1) 



of a SMBH will be tidally disrupted. In this equation A/bh 
is the black hole mass, m, and r» are the stellar mass and 
radius, and rj is a stellar structure constant which is of or- 
der unity for main sequence stars. Detailed calculations find 



77 = 0.844 for n = 3 polytropes, for example (Diener et al. 
1995). During the disruption event, roughly one half of the 
stellar mass is ejected from the system, while the other half 
remains bound (Rees 1988) but with the change in gravita- 
tional potential across the star producing a wide spread of 
specific orbital energies. 



GMbh 



(2) 



Here rp is the pericenter distance of the star on its orbit 
around the black hole. The bound stellar matter quickly 
expands to the point where hydrodynamic forces can be ne- 
glected and the gas follows roughly Keplerian trajectories, 
but upon return to pericenter the gas streams dissipate en- 
ergy in shocks and form an accreting torus. The character- 
istic mass return rate is 



Mr 



1 m* 

3ir 



-5/3 



(3) 



(Phinney 1989), with the return time for the most tightly 
bound debris being 



tr 



2-R 



3/2 



\GMbh 



1/2 



(4) 



(SQ09). This simple dynamical picture has been confirmed 
as largely accurate by numerical hydrodynamic simulations 
(Evans & Kochanek 1989). The radiative properties of the 
flare have been modeled by Loeb & Ulmer (1997), and most 
recently by SQ09, whose model finds that mass infall rates 
can exceed the Eddington limit for days to months after dis- 
ruption. This super-Eddington phase may produce an out- 
flow with supernova-like optical luminosities. After super- 
Eddington infall ceases, emission is dominated by a thermal 
component from the accretion disk, which also photoionizes 
the unbound stellax debris, producing broad emission lines. 
Luminosity from the disk fades with a decay time of order 
months to years, and is peaked in the UV and soft X-ray. 
Upcoming transient surveys will dramatically expand the 
number of TDEs available for study, with the Pan-STARRS 
37r survey expected to flnd ~ 10 per year, rising to ~ 200 
if the SQ09 predictions about super-Eddington outflows are 
correct. LSST is expected to find ~ 100 per year, or up to 
~ 6000 if the hypothesized super-Eddington outfiows can be 
seen. In this paper we consider a black hole mass range from 
W'^Mq to 10* M0. For Mbh > 10* Mq the Schwarzschild 
radius of a black hole exceeds the tidal disruption radius of 
main sequence stars, meaning that stars will bo swallowed 
whole rather than disrupted by non-spinning SMBHs. How- 
ever, for a Kerr black hole rt is angle-dependent, and for high 
values of spin, black holes as large as ~ 7 x 10* Mq can still 
tidally disrupt main sequence stars which approach from an- 
gles near the equatorial plane (Boloborodov ct al. 1992). For 
simplicity we consider SMBHs up to but not above 10* Mq. 

In this paper we follow the prescriptions of SQ09 and 
Strubbe & Quataert 2010 for disk and super-Eddington out- 
fiow luminosity. Specifically, we model the super-Eddington 
outflow as a sphere with photosphere radius 



Rph ~ 10/out/v ^ 



(5) 
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photosphere temperature 



Tph ~2 X 10 



fou 



-5/12 



(6) 



and time of peak emission equal to the larger of ti and the 
time when the outer edge of the photosphere becomes opti- 
cally thin, 

W - 1 f!i^f.-'^'Ml^'Rl%,rnrr-'^' days. (7) 
Here the outflowing wind velocity is given by 



Vv>, = /v 



(8) 



with /v a free parameter of fiducial value 1, and the total 
mass flux in the outflow is 



A'/out — font Ml, 



(9) 



with /out a parameter taken to be 0.1. MEdd is the SMBH's 
Eddington-limited mass accretion rate assuming an ac- 
cretion efficiency of 0.1, Me = Mbh/(10''Mq), Rs = 
2GMBH/c^ and Rp,aRs = rp/{3Rs). 

Again following SQ09, we approximate the disk as a 
geometrically slim multicolor blackbody extending from the 
innermost stable circular orbit, -Risco, to 2r-p. Its tempera- 
ture profile is given by 



„4 SGMbuMJ^ 



(10) 



1 /l 3 / lOM,it:s 

2 V4 2HMEddi? 



with / = 1 — \/ Risco /R, and -Risco a function of black hole 
spin. 

The existence of suitable emission or absorption lines 
for kinematic identification of a recoiling TDE is an open 
question. There is widespread agreement in the literature 
on the existence of emission lines in the immediately un- 
bound, photoionized stellar debris (Bogdanovic et al. 2004, 
Strubbe & Quataert 2009, Strubbe & Quataert 2010). Some 
of these spectral features are potentially useful for distin- 
guishing TDEs from supernovae and other transients (Kasen 
& Ramirez-Ruiz 2010). Unfortunately, the velocity spread in 
this debris. At; ~ {2Ae)^^^ , can easily reach a large fraction 
of the speed of light, and details of the Doppler broaden- 
ing will be determined by the inclination of the star's ini- 
tial orbit around the black hole, among other unknown pa- 
rameters. A more promising candidate seems to be absorp- 
tion lines formed when outflowing material processes con- 
tinuum radiation from the accretion disk around the black 
hole (Strubbe & Quataert 2010). Most of these lines are in 
the UV part of the spectrum, although weaker hydrogen or 
helium lines may exist in optical bands, and the presence of 
a soft X-ray, power law tail (as has been observed in tidal 
flares detected by ROSAT and GALEX - see Komossa 2002, 
Gezari et al. 2008) could produce absorption lines between 
1 and 10 keV as well (Strubbe & Quataert 2010). If the 
wind launching speed does not vary too much in time, these 
lines will be sufficiently narrow (with thermal broadening 
~ 30 km s~^) to make the velocity difference between the 



outflowing wind's photosphere and the host galaxy measur- 
able. However, the bulk outflow velocity may itself be quite 
large, in which case determination of the black hole veloc- 
ity would require extremely precise wind modeling. Alter- 
natively, if /v ^ 1, a recoiling black hole's velocity could 
dominate the kinematic offset between the photosphere and 
the galaxy, and make wind launching speed variations small. 
A limiting case of this regime is a different model for the 
super-Eddington phase of mass return (Loeb & Ulmer 1997), 
in which radiation pressure from the disk isotropizes return- 
ing gas streams and supports a quasi-spherical cloud of dis- 
rupted matter. In this model, measurement of the absorption 
offset from host galaxy lines would precisely identify which 
TDEs were caused by recoiling black holes. Although the 
validity of either of these models is at the moment unclear, 
recent observations of TDEs found in SDSS data (van Velzen 
et al. 2010) may be more compatible with the predictions 
of Loeb & Ulmer 1997. Further validation of this model (or 
a low /v version of SQ09) would indicate the feasibility of 
kinematic identification of recoil-induced TDEs. 

For a recoiling SMBH, there are two sources of stars 
to tidally disrupt; unbound stars encountered in its passage 
through the host galaxy, and the small cloud of stars in 
orbits that remained bound during the recoil event. Both 
of these sources are considered in the next subsections. For 
simplicity, in the remainder of this paper we set rj = 1, m* = 
Mq , and r* — Rq (a conservative assumption - see Tremaine 
& Magorrian 1999 for how generalizing to a more realistic 
stellar mass function can increase tidal disruption rates by 
a factor ~ 2). 



2.3 Host Galaxy Structure 

Like many of the stationary SMBH tidal disruption rate 
papers (Ulmer & Syer 1998, Tremaine & Magorrian 1999, 
Wang & Merritt 2003) we consider the Nuker surface bright- 
ness parametrization (Lauer et al. 1995) 



I{r) = /b2 



B-r / r 
rh 



i + (^ 



(11) 



which was originally developed to model the surface bright- 
nesses of nearby ellipticals and bulges resolved at the parsec 
level with the Hubble Space Telescope. In Equation (|ll|l , rb 
is the break radius at which the shallow inner power law of 
I{r) oc becomes the outer power law, /(r) oc r"^. The 
strength of the break is determined by the dimensionless 
coefficient a, and Ih is the surface brightness at the break 
radius. The profile can be Abel inverted to yield a volume 
density profile with power law exponents of 7 ~ F + 1 and 
P ^ B + 1 when r <^ and r 3> ^b, respectively. The 
Nuker parametrization has many well-known surface bright- 
ness models as special cases (Byun et al. 1996). 

The Nuker galaxies are roughly divided into two cate- 
gories based on the steepness of the surface brightness profile 
in the centermost regions: core galaxies and cusp galaxies. 
It has been suggested (Merritt & Milosavljevic 2005) that 
core galaxies are formed when the inspiral and merger of 
two SMBHs ejects stars from the host galaxy's center via 
3-body interactions, scouring a core from the stellar profile 
- though subsequent star formation can rebuild the central 
parsecs into a cusp in the case of gas-rich (so-called wet) 
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mergers (Kormendy at al. 2009, Blecha et al. 2010). For this 
reason, we will take core galaxies as those which best rep- 
resent the immediate post-recoil state of the stellar profile 
in gas-poor (so-called dry) mergers. We will consider cusp 
galaxies also, as they arc likely to be relevant in wet merger 
scenarios. Most core galaxies in the Nuker sample tend to 
possess large central black hole masses Mbh 10* M©, as 
determined by the Mbh — cr relation). To examine a wider 
range of galaxy masses, we do not directly use observed sam- 
ples of galaxies but rather create a simple mock catalog of 
galaxies in which we simulate the black hole trajectories. 
Important scaling relations which we require for all galaxies 
in our mock catalog are the Mbh — cr* relation (Tremaine et 
al. 2002), 



cr* = 220km s 



Mb 



2 X 108 Mp 



(12) 



and the bulge-black hole mass relation (Marconi & Hunt 
2003) 



M* = 4.06 X 10'" Mq ( 



V 108 M© 



(13) 



with M, being total stellar mass and cr, the bulge velocity 
dispersion. We also fit our galaxies to empirically determined 
calibrations of SMBH influence radii (Merritt et al. 2009): 



r-infl = 35 pc 



nnfl = 22 pc 



/ Mbh 

V 108 Mq 



Mb 



(l08M© 



(14) 



(15) 



The first of these relations holds for core galaxies and the 

second for cusp galaxies. Here the influence radius, rinf is 
the radius within which there is a mass in stars equal to 
2Mbh. 

The most relevant aspect of the Nuker parametrizar 
tion for our purposes is the broken power-law in the stel- 
lar density profile; as we shall see, the trajectories of kicked 
black holes are strongly affected by central density gradi- 
ents. While the most self-consistent way to proceed would 
be to numerically deproject scaled brightness profiles 7(r) 
into luminosity density profiles j(r), we suspect this would 
fail to capture an important part of the physics of black hole 
recoil. Initially radial orbits in spherically symmetric poten- 
tials make multiple passes through galactic center, whereas 
initially radial orbits in nonspherical potentials will gener- 
ally receive torques that increase their pericenter distance. 
This is of importance to this paper because orbital energy 
loss for recoiling black holes in dry mergers is dominated by 
stellar dynamical friction during passes through the densest 
central regions (Gualandris & Merritt 2008). Using a spher- 
ical potential would produce purely radial orbits, artificially 
shortening the wandering life time of the black hole, and 
decreasing the number of offset TDEs it can produce. To 
avoid this problem, we use an ad hoc density-potential pair 
that is designed to approximate the most important dynam- 
ics of a Nuker profile but which is also easily generalizable 
to nonspherical geometries. Specifically, we create a set of 
spheroidal isodensity surfaces 



p»(m) = < 



Kim ' m < mb 
K2m~^ rub < mmax 
m ^ mmax 



using the axisymmetric pseudoradius 



2 r z 
m = ^ + ^. 



(16) 



(17) 



In these equations r and z are standard cylindrical coordi- 
nates, and the isodensity spheroids have dimensionless axis 
ratios a and c (a > c). The so far free parameters 7, /3, 
and mb are calibrated using the corresponding values in 
the spherical spatial deprojections of the most recent Nuker 
sample (Lauer et al. 2005). None of those three parame- 
ters vary strongly with A/, once the core/cusp dichotomy 
is accounted for, so we treat them as having Gaussian dis- 
tributions about their measured means, with variance also 
calibrated off the Lauer et al. 2005 sample. Ki and K2 arc 
chosen so that the density profile is continuous at the break 
pseudoradius, mb, and mmax is selected in combination with 
Ki,K2 to both normalize the total stellar mass, M*, and to 
reproduce the correct value of rinfl. The final free param- 
eter of our mock galaxy catalog is ellipticity e — 1 — c/a. 
Fortunately for our purposes, Lauer et al. (2005) measured 
isophote ellipticitics for the galaxies in their sample, which 
we use to sample ellipticity parameter space (again, due to 
its lack of variation with M* , we treat ellipticity as Gaussian- 
distributed) . 

The distinction between cores and cusps deserves more 
consideration, however. Because SMBHs above 10* Mg con- 
sume stars whole rather than tidally disrupting them, the 
majority of galaxies of interest to us fall into the meiss range 
represented primarily by cusps. Although scaling relations 
for core galaxy parameters exist (Faber et al. 1997), it is 
not clear how reliably they can be extrapolated an order of 
magnitude below the smallest core galaxies in the samples 
which they are based on. For our mock catalog, we consider 
the mergers of initially cuspy galaxies, but treat separately 
two different limiting scenarios: 

• In the first, the tendency of a binary SMBH inspiral 
to scour a core is outmatched by star formation, and a nu- 
clear cusp is preserved. In this pure cusp scenario, all free 
parameters {7, /3, rb, e} are calibrated off the sample of cusp 
galaxies. We assume in this case that a fraction /g of the 
baryonic mass of the bulge mass is in the form of a gas disk, 
which we describe in more detail in §2.4. 

• In the second, we consider a SMBH binary which is 
successful at scouring a core, either because its progenitor 
merger was dry or because free gas was consumed or expelled 
prior to the late phase of the hard binary's evolution. For 
this case we calibrate the constant value of /3 and the initial 
values 70, rb,o off the cusp sample, but calibrate e and final 75 
off the core sample, and manually "excavate" a mass deficit 
A A/ — 2Mbh to determine the final break radius rb,f. The 
final value for the new, cored break radius is 



rb,f = 



47rKi(l-e2)l/2 



.3-^ 

'^b.O (3-^)(3-7o) 



-Tn 



1/(3-/9) 



(18) 



'b,0 (3-7f)(3-/3)) 

In practice, this tends to increase the break radius by a 
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factor of a few. Here we have used the elhpsoidal eccentricity 
e = \/l — [c/a)'^. In this second scenario, we set the remnant 
gas fraction /g = 0. 

With our galaxy models fully determined, we can then 
integrate our density profiles numerically to calculate other 
relevant quantities such as the stellar potential and forces. 
Here we use the standard method of homoeoids (Binney & 
Tremaine 2008, Section 2.5). Due to the difficulty involved 
in even numerical calculation of a two-integral distribution 
function (Hunter & Qian 1993), particularly when an analyt- 
ical, closed-form potential is lacking (as is the case hero) we 
are forced to compute the velocity dispersion in the spherical 
limit and generalize by substituting m for r; this introduces 
modest inaccuracy into the dynamical friction and gravitar 
tional focusing formulae used in §2.4. The spherical limit of 
this broken power-law profile has a velocity dispersion 



a{r) = 2nGK2i 



2r 



3-/3 



r(3-7)(l + /J) 



(19) 



-I- 



for r > rb, and 



2r; 



,3-/3 



(3-/3)(/3-l) r(3-/3)(l+/3) 



fT^(r) = 2nGKir 



(3-7)(l + /3) 
2 



(3-/3)(/?-l) (3-/?)(l+/9) 



(20) 



+ 



27rGj^ir^-T 
(3-7)(l-7)' (3-7)(l-7) 



for r ^ Tb. 

The stellar bulge population dominates gravitational ef- 
fects for low velocity kicks, but higher velocity kicks carry 
SMBHs into regions where the dark matter halo potential 
becomes important. We model the density profile of dark 
matter using a standard NFW profile (Navarro et al. 1997), 
given by 

Po 



Pnpw ~ 



{r/a,){l + r/a,y 



(21) 



Here the scale distance Os and density factor po are de- 
termined by assuming a concentration of 10 and truncat- 
ing the NFW profile at r2oo, the virial radius at which 
Pnfw = 200pc, with pc the cosmological critical density. 
We also normalize the total dark matter mass Mtot using 
the latest calibration of the MBH-Mtot relation (Bandara et 
al. 2009), 



Mbh = 1.51 X 10** 



Mtc 



V 1013 Mo 



(22) 



Our assumption of spherical symmetry for the dark matter 
halo should be considered conservative, since nonspherical 
potentials extend the wandering time of the SMBH. 

Our general strategy for galaxy modeling is to select a 
fiducial value of SMBH mass, a parameter set {7, 13, rb, e} in- 
formed by the Laucr et al. (2005) sample, and then to use the 
scaling relations described in this section to self-consistently 
find other parameters so that integrated quantities like po- 
tential or velocity dispersion can be numerically computed. 



The fiducial values of Mbh wc use are 10^ 10** ^ 10^ 10^ ■^ 
and 10* M0. To parallel the SQ09 calculation we also use the 
black hole mass function inferred by Hopkins et al. (2007b). 

In minor mergers, the dynamical friction timescale for 
satellite infall is expected to be greater than the Hubble 
time (Wetzel & White 2009), meaning that we only need to 
consider the mass range of g = 0.1 — 1. Simulated merger 
rates for this mass ratio range in galaxies with stellar mass 
M* > Afmin have been matched (Hopkins et al. 2009) to the 
analytic fit 



dt 



=0.04 1 + 



Mo 



(23) 



x(l + ^)«"— ' Gyr-\ 
with z the redshift of the galaxies, Mo = 2 x 10^" M©, and 



/3(M„,i„) = 1.5 -0.25 log 



Mo 



(24) 



Wo use this merger rate in combination with the volumet- 
ric SMBH mass function and the results of our orbit in- 
tegrations (described below) to compute the total rate of 
events observable by LSST. For all calculations in this paper 
we assume a standard ACDM cosmology with Qa = 0.73, 
= 0.27, and Hq = 71 km s"i Mpc"i (Spergel et al. 
2007). 

2.4 Interactions with the Galaxy 

To model the disruption rate of unbound stars, we evolve 
the SMBH's trajectory under the influence of gravity and 
dynamical friction through galaxies with the axisymmetric 
Nuker density profiles described above. Wo use the fifth- 
order Dormand-Prince method (with an embedded fourth- 
order Runge-Kutta integrator for adaptive timestepping) de- 
scribed in Numerical Recipes, Chapter 17. The effects of dy- 
namical friction are approximated with the Chandrasekhar 
formula (Chandrasekhar 1943), 



Fdf : 



-7(M) 



AnpiGAlBnf VBH 



(25) 



with cr the local velocity dispersion and p the local density 
of the medium causing the drag. For a collisionless medium, 
such as a stellar population. 



7dry(M) 



ln(A) 

M2 



erf 



M_ 
V2 



Me 



(26) 



where M = vbh/ct is the Mach number. The Coulomb log- 
arithm can be fit numerically (Escala et al. 2004), and for 
the case of SMBHs on radial orbits, is well fit by a value of 
InA = 2.5 (Gualandris & Merritt 2008). Dynamical friction 
is the force which ultimately causes the kicked black hole to 
settle back into a near-stationary position in the center of 
its host galaxy, on timcscalcs ranging from 10" to lO** years 
(Madau & Quataert 2004, Blccha & Loeb 2008). The Chan- 
drasekhar formula is derived assuming a uniform and infinite 
background of stars, and it is not immediately clear how ap- 
propriate that is for a steep density profile in galactic nuclei, 
or for a black hole massive enough to excite a response in the 
stellar population. The applicability of the Chandrasekhar 
formula to bound, recoiling black holes has been considered 
before (Gualandris & Merritt, 2008), and for appropriately 
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chosen In A it was found to be fairly accurate until the point 
when the mass interior to the black hole's apogalacticon is 
of order Mbh. After this, coherent oscillations develop in 
the stars interior to the black hole's trajectory, and dynam- 
ical friction is found to become dramatically less effective at 
removing the black hole's orbital energy. We terminate our 
calculations at the onset of this orbital phase, both because 
our trajectory approximation would become quite inaccu- 
rate and also because TDEs caused by a slow-moving SMBH 
near the center of a galaxy would not be distinguishable from 
those caused by a stationary black hole. During these calcu- 
lations, we neglect the extra "core scouring" caused by black 
hole recoil (Gualandris & Merritt 2008) . The stellar popula- 
tion in the galactic center responds to a moving SMBH by 
expanding, with a mass of stars equal to a few A/bh being 
displaced from the galactic center for kicks close to escape 
velocity (and the effect is reduced for slower ones). Neglect 
of this effect likely reduces the SMBH wandering time and 
causes us to underestimate the total number of TDEs per 
galaxy merger, but probably not by much, as axisymmetry 
of the stellar potentials prevents the SMBHs from returning 
exactly to the center of their host galaxies where core scour- 
ing is most relevant. We highlight that dynamical friction re- 
moves the most orbital energy during passages through the 
densest regions of the SMBH's trajectory. Therefore trajec- 
tories with nonzero angular momentum (due to an axisym- 
metric potential) last longer before settling back into the 
galactic center than would center-crossing ones in spherical 
geometries. 

In a dry merger it is sufficient to consider dynamical 
friction off stars and not gas. This regime could also apply 
to wet mergers where the gas is used up in star formation 
(while the SMBH binary is stalled) or dispersed in binary 
quasar feedback. We identify both of these scenarios with 
our "excavated core" galaxies. If significant quantities of 
gas survive until the recoil phase of the merger, however, 
it is necessary to consider the effects of gas dynamical fric- 
tion, which would apply more to our "pure cusps." Previous 
work (Blecha et al. 2010) has indicated the effect of left- 
over gas is to decrease black hole wandering times, reducing 
the observable number of offset TDEs. To quantify this ef- 
fect, Equation (|26p still applies; we simply need to use gas 
density rather than stellar density for p, substitute a local 
sound speed Cs for a, and modify the dimensionless param- 
eter I{M) (for gas, M — ubh/cs). The new dimensionless 
functions are 



-^subsonic (A-^) — - In 



in the subsonic regime, and 



1 + M 

1 - M 



/supersonic (M) = ^ ( ^ ^ 



In 



M 



VBut 



(27) 



(28) 



in the supersonic regime (Ostriker 1999). However, these 
formulae have been shown to overestimate gas dynamical 
friction in the slightly supersonic regime, so we adopt the 
prescription of Escala et al. 2004 and use the Chandrasekhar 
formula for I{M), with InA = 2.5 for M ^ 4.7 and InA = 1.5 
for M < 0.8. We follow the prescriptions of Blecha & Loeb 
(2008) and assume that most of the gas in the galaxy has 
settled into a disk, which we align with the oblate plane of 
the galaxy. We employ a slightly less complicated version of 



their model, as only two of their four disk zones are relevant 
for our dynamical modeling: zones III and IV (zones I and 
II only exist in the presence of a central SMBH). Zone HI, 
the portion of the disk influenced by the SMBH potential 
before the recoil kick, is truncated on its inner edge at the 
kick radius, 



GMbk 



(29) 



and transitions to zone IV at r = rinfj. Zone IV is an ex- 
ponential disk with scale rdisk- The disk surface density in 
zone HI will be 

f2 X 1/3 



(30) 



(31) 



while in zone IV, the surface density is 

_ A:fdisk(> Hnfl) e'''^'' 

^ 27rrdisk(nnfl +rdisk)e-'-i"ti/'-disk 

Here we take the viscosity parameter a = 0.1, and set 
the Toomre parameter Q = 1 (Toomre 1964) under the 
assumption that star formation feedback roughly balances 
cooling, leaving the disk marginally stable. The scale dis- 
tance rdisk is found by requiring continuity between zones 
HI and IV: Eiii(rinfl) = Eiv(nnfl). The accretion rate Ala 
can be found by mass normalization of equation (|30p so that 
Mm = 2/gasMBH: 

Ma = ( f^^M^y aG'^\ (32) 
V4(rinfl-rk)/ 

Densities in both disk zones decay exponentially with height 
z, with scale height 

(33) 



mil = —r 



taken from Blecha & Loeb 2008. In zone IV, the scale height 



MaKn 



(34) 



is solved for using the identities Ma ~ SnaCshT, and Cs = 
(GMa/ay^'-'- Here Kn is the epicyclic frequency, and is cal- 
culated from the numerically integrated potentials of the 
isodensity shells in Equation (|16|l . 

In our models we consider values of /g (gas as a fraction 
of total baryonic mass) of and 0.3. The latter value is taken 
as a conservative upper limit for remnant gas fraction at the 
time of black hole merger, as self-consistent hydrodynamical 
simulations (Mihos & Hernquist 1996) have shown that ^ 
50% of the initial gas fraction, /g^ in a merger is expelled or 
converted into stars by the time of black hole coalescence. 
Observation indicates that /gj < 0.6 for Mbh > 10® M© 
at low redshift (Hopkins et al. 2009, Figure 7), so /g — 
0.3 is a conservative case, likely to result in SMBH orbits 
which decay somewhat faster and produce fewer offset tidal 
disruptions than in more general wet mergers with smaller 
/g- 

At each point along the SMBH's trajectory we consider 
an instantaneous "tidal disruption cylinder" of length UBuAf 
and radius equal to the gravitationally focused tidal disrup- 
tion radius. This lets us simply calculate instantaneous tidal 
disruption rates along the trajectory, 

2GMbh \ 



A^u = ptV-rrrt 1 -I- 



rtv^ 



(35) 
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which can be integrated to get a time-averaged TDE rate, or 
A^'tde, the total number of stars disrupted per recoil event. 
Here v = \A'^^7+"^' with a given by Equations (|19|l and 



2.5 Interactions with the Bound Cloud 

The initial size of the bound cloud is determined by the 
magnitude of the received kick, and can be approximated 
as encompassing all stars within rk. The mass of the bound 
cloud is found by KM08 to be a fraction fb of the black hole 
mass, where 

3-7 

(36) 



Here 7 is the same as in Equation (|16p . rinfl is the influ- 
ence radius, the interior of which contains a mass in stars 
twice Mbh, and -^(7) = 11. 67"^'^^. For most cloud sizes the 
disruption rate of bound stars will be determined by reso- 
nant relaxation into the SMBH's empty loss cone, exponen- 
tially depleting the population of stars inside on a timescale 
r « 3.6GM|H/(«k?n,) (KM08). In practice, the e-folding 
time is at least an order of magnitude below 10^" years for 
most of the 10^ Mq and 10^'"' Mq black holes which escape 
from their host galaxies, strongly suppressing the averaged 
intergalactic TDE rate. 

One exception to this picture is if rk ~ Hnfl; in this case 
nonresonant relaxation could become important, and res- 
onant relaxation alone will significantly underestimate the 
TDE rate. This regime is of minimal significance for this 
paper, however, since small kicks are likely to produce few 
spectrally and no spatially offset flares. A more significant 
exception is for relatively low-mass SMBHs, which can reach 
an energy relaxation timescale in less than their wandering 
time. Relaxation will eventually allow the cloud to expand 
in radius (O'Leary & Loeb 2009), changing the time evolu- 
tion of the tidal disruption rate from exponential depletion 
to oc t"^/^ (O'Leary & Loeb 2011). Therefore, we adopt 
KMOS's prescription for resonant relaxation, 

iV,^CR,R,(7)^^^^^/be-*/^ (37) 
ln(rk/rt) T\ 

when t < tr, but transition to iVb oc t~^^'^ at later times. 
The energy relaxation timescale tr is taken to be (O'Leary 
& Loeb 2011) 

5/4 / ^ \ 1/4 



tr 



W yrs 



Mbh 
105 Mp 



rk 

rinfl 



(38) 



This power law disruption rate is only relevant for Mbh < 
10*' '^ Mq, but for lower mass SMBHs we transition to power 
law depletion after an energy relaxation time. For both sce- 
narios, the initial disruption rate is 

-2 



Nb 



1.5x10" 



Mbh 



10 pc 

-1 

yr" 



(39) 



Vl03 km s-i. 

One uncertainty is the resonant relaxation coefficient Crr, 
found by KM08 to have a value of 0.14 for 7=1. Since the 
spatial power law exponents for core galaxies are close to 
1, we adopt this value, though it is less well motivated for 
cuspier galaxies. 



We also consider growth of the bound cloud by cap- 
ture of members of binary star systems. This three-body 
interaction is treated in the same way as tidal disruption of 
unbound stars, except instead of a stellar tidal disruption 
radius we use an "orbital tidal disruption radius" , given by 



rt,o = flbin 



V 2m. 



1/3 



(40) 



where abin is the binary semimajor axis. While one member 
of the binary is ejected at high velocities (Hills 1988), the 
other is bound to the black hole, with apoapsis rmax given 
by 



GMbh 



(Yu & Tremaine 2003), with Wcjoct ~ 145km s"^ (Hills 1988). 
To calculate the rate of these captures, we assume Opik's 
Law (Opik 1924), a fiat distribution of binary semimajor 
axes a in units of log(a), between amin and Umax- Following 
Vereshchagin et al. (1988) and Kouwenhoven et al. (2007) 
we adopt Umin — 5Rq and amax — 5x 10^ -Rq. We only con- 
sider captures with rmax < rt,c, with the cloud's tidal radius 
conservatively given by rt^c ~ rinfl. This refill mechanism is 
in principle capable of counterbalancing losses due to tidal 
disruption and evaporation from the cloud. Without a refill 
source, resonant relaxation into the loss cone will normally 
cause the population of the bound cloud to evolve due to 
A^TDE oc N, leading to a population (and TDE rate) de- 
pleted exponentially in time. We can roughly see the effect 
of stellar capture into the bound cloud if we assume 



TV 



~kN + m, 



(42) 



with k the average frequency with which bound stars evap- 
orate or are scattered into the loss cone and m a time- 
averaged capture rate. This differential equation has the so- 
lution 



iV(t) = 7V(0)e-** -1-^(1 
k 



(43) 



By itself, the resonant relaxation process will deplete the 
bound cloud, but 3-body capture allows the number of stars 
in the cloud to asymptotically approach a nonzero value. If 
the time-averaged binary capture rate is high enough (i.e. 
if m/k > N{Q)) the size of the cloud would even grow over 
time. The importance of this effect is determined for each 
galaxy/kick velocity pair. 

A final consideration is stability of the bound cloud 
to perturbations. Analytically, it seems unlikely that in- 
teractions with unbound stars will eject significant num- 
bers of bound stars from the cloud: if the cloud stars are 
bound to the black hole with typical energy -Bbind ~ —v^rrit,, 
and during encounters with unbound stars a change in en- 
ergy A_E ~ Gml/rp is available (where rp is the clos- 
est approach of the two stars), encounters must be within 
7" < rp ~ rk(mt/MBH). For a IO'^Mq black hole on typical 
trajectories, this works out to at most ~ 1 unbound stars 
making close enough approaches to eject a bound star during 
the SMBH's passages through the bulge. 
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2.6 Observability of Recoil-Induced TDEs 

To translate the total recoil-induced TDE rate into a rate 
of identifiably recoil-induced TDEs, it is necessary to con- 
sider observational constraints. LSST's rapid cadence, high 
sensitivity and thorough sky coverage make it an ideal sur- 
vey to detect disruption flares - as mentioned in §2.2, LSST 
could detect up to thousands of TDEs per year. LSST's lim- 
iting g-band magnitude is 25 (LSST Science Collaborations 
2009); because of LSST's short cadence we assume any flares 
brighter than that will be detected. The detectability of a 
spatial offset will depend on how well the TDE centroid 
can be distinguished from the host gala^cy centroid after 
photometric frame subtraction. For LSST the expected dif- 
ferential astrometric precision will be ~ 0.7"/SNR (LSST 
Science Collaborations 2009). Using the LSST Science Man- 
ual's prescription for SNR"-*^ — atot = \/ o-|ys + ""rand ' '^^ 
infer astrometric precision by calculating the signal to noise 
ratio for each event in our sample. We also calculate the 
rate of spectroscopically identiflable flares associated with a 
recoiling SMBH. Although UV spectroscopy would be ideal, 
soft X-ray spectrometers - SXS, for example, on the planned 
ASTRO-H mission (Takahashi et al. 2010), expected to be 
operating contemporarily with LSST - should be able to 
identify the absorption lines discussed in §2.2, if they exist 
with suflicient equivalent width. To investigate this possibil- 
ity, we consider a fiducial case of absorption lines at 10 keV, 
observed by SXS foUowup with an energy resolution of 7 
eV. If the outflowing wind can be accurately modeled, these 
lines would allow black hole velocities down to ~ 200 km s~^ 
to be spectrally resolvable. As mentioned before, it is not 
clear that the super-Eddington phase of accretion will pro- 
duce winds in which a ~ 200 km s~^ offset is detectable, so 
predictions of kinematic offsets should be regarded as some- 
what hypothetical. Because spatial and kinematic offsets are 
angle- dependent, we average the observable quantities over 
all inclination angles for the host galaxy. 



3 TDE RATE 

Using the potentials and frictional forces described above, 
we integrate the trajectories of flve different black hole 
masses at nine different kick velocities and seven inclination 
angles in galaxies with eighteen different possible permuta- 
tions of mass-independent structural parameters, for a total 
of 11340 runs (the flnal factor of two comes from wet vs dry 
mergers). During preliminary test runs, a very weak depen- 
dence of the wandering time on /3 and 7 was apparent (once 
variation in /3 and 7 due to the core/cusp dichotomy is al- 
lowed for), so we set those quantities equal to their average 
values. Among the remaining structural parameters, we only 
varied e and rh,o- 

We terminate our trajectory calculations after a Hubble 
time, if the black hole has left the stellar bulge (and its 
attendant sources of friction) with escape velocity, or upon 
the onset of the "Phase H" orbital oscillations of Gualandris 
& Merritt (2008), discussed in §2.4. 

To calculate the total observable rate of TDEs due to 
recoiling black holes, F, we use a modifled version of Equa- 



tion (31) in SQ09. Speciflcally, 

d7(r, rp) 
X M rdlnrp 
dlnrp 

Here dmax('"p) represents the maximum comoving distance 
a TDE flare with given pericenter rp can be seen at, using 
a 25 AB g-band magnitude limit. LSST will detect flares 
at cosmological distances, so it is necessary to employ a K- 
correction, which has a modest impact on dmax. The rate 
7(r, rp) is integrated over inclination/azimuth angles and 
galaxy properties, and is the rate at which either TDE flares 
are produced at a distance r by TDEs with rp. Further cuts 
are added to the integrand to calculate the rate at which 
observably spatially offset TDE flares, or observably kine- 
matically offset TDE flares are produced, using the crite- 
ria described in §2.6 (with an average over azimuthal an- 
gles to account for projection effects). In our average over 
galaxy properties, we give MsH-dependent weights to the 
"pure cusp" and "excavated core" scenarios from §2.3. These 
weights are determined by binning the the Lauer et al. 2005 
galaxy sample and computing the fraction of cusps and cores 
in each mass bin (with the small minority of intermediate 
cases taken as 50% core, 50% cusp). 

To calculate dmax for super-Eddington flares we use 
Equations ([5)l, (|6]), and ((Tjl, while to do the same for disk 
emission we use Equation (llOp . For simplicity, we neglect the 
less important emission from photoionized, unbound disrup- 
tion debris, noting that this is a conservative approximation. 
For disruptions from the bound cloud, resonant relaxation 
slowly diffuses stars across the loss cone in phase space, 
meaning that nearly all bound TDEs will have rp ~ rt. Un- 
bound stars will have a wider variety of rp, but the geometry 
of gravitational focusing will bias them towards rp « rt as 
well. For these reasons we simplify Equation (|44p by tak- 
ing rp = rt . This approximation produces slightly more disk 
emission (due to physically larger disks), and signiflcantly 
less luminous super-Eddington flares, than the flat distribu- 
tion of TDEs across Inrp assumed in SQ09. Consequently, 
our results show a much less pronounced difference in the 
observable TDE rate between the disk emission and super- 
Eddington outflows cases. 

We then interpolate the results of these trajectory cal- 
culations over three different black hole physics scenarios, as 
discussed in §2.1. In the first scenario, a lack of free gas dur- 
ing the SMBH inspiral leaves the spin vectors of the SMBH 
binary randomly aligned with each other, producing a top- 
heavy kick distribution and a high average value of dimen- 
sionless spin (a = 0.73). The other two scenarios involve wet 
mergers with warmer and cooler gas, producing spin vectors 
aligned to within 30 and 10 degrees, respectively, and rem- 
nant mean spins of a = 0.88 and a = 0.90. Because the 
disruption of stars by 10* M© black holes is so sensitive to a 
(a > 0.92 required), we bracket these fiducial assumptions 
about remnant spin with a — and a = 0.95 cases. We also 
consider two different cases of tidal disruption physics; one in 
which the super-Eddington mass outflows proposed by SQ09 
exist (for simplicity we take their canonical case of /v = 1 
and /out — 0.1), and the other in which they do not. In the 
latter, optical emission is limited to the Rayleigh- Jeans tail 
of the newly-formed accretion disk. 
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Figure 1. The importance of nonsphericity on SMBH trajecto- 
ries. All the lines represent trajectories of 10^ Mq SMBHs kicked 
at 400 km s~^ in a core galaxy. The green solid line is the course 
of a SMBH kicked in the plane of the galaxy's two semimajor 
axes, while the blue dot-dashed line represents a SMBH kicked 
15 degrees above the plane, and the black dotted line a SMBH 
kicked 45 degrees above the plane. 

4 RESULTS AND DISCUSSION 

Figures [1] and [2] illustrate the importance of nonspherical 
potentials for the lifetimes of wandering black holes. In the 
axisymmetric stellar potential that we employ, stars initially 
on radial orbits will quickly acquire angular momentum un- 
less they lie on a principal axis of the stellar ellipsoid or in 
the equatorial plane. The latter is true of the 0°, green orbit, 
which is seen in Figure 1 to decay somewhat faster than the 
blue, 15° orbit and dramatically faster than the 45°, black 
orbit. For the dry mergers illustrated here, the variation in 
decay time is due entirely to differences in stellar dynamical 
friction, which is the strongest at orbital pericenter. In the 
wet merger scenarios we considered, recoils in the plane of 
the gas disk are very quickly damped out, but axisymme- 
try in the stellar potential still affects orbital lifetimes for 
other inclination angles. Figure[2]illustrates the torques that 
act on orbits out of the equatorial plane, and in the inset 
we can see that those torqued orbits are able to avoid close 
pericenter passages, explaining their longevity. 

The large density variations across the SMBH host 
galaxies can be seen in the unbound stellar disruption rate, 
displayed in Figure [3l which spans many orders of magni- 
tude. The spikes are located at passages through the galac- 
tic center, and their increasing magnitude with each cycle 
arises from the inverse relation between gravitationally fo- 
cused cross section and velocity. The thick lines, representing 
the disruption rate of bound stars, are much more constant 
in time, although for lower-mass SMBHs these become ob- 
servably depleted, as discussed earlier. In Figured we plot 
the unbound stellar disruption rate versus radial distance, 
and can clearly see the break in the stellar density profile 



Figure 2. The same black holes as in Figure 1, now viewed in two 
dimensions. The inset plot zooms in on the central 100 parsecs 
to highlight the lack of center crossings for SMBHs ejected at 
nonzero inclination angles. 

near 7 pc. Scatter in Figure |4] is due to the combination 
of different velocities and different inclination angles during 
pericenter passage. The binary capture effect hypothesized 
in §2.5 was seen at low levels but found not to contribute 
substantially to bound cloud sizes or disruption rates. 

The total number of disruptions per merger, A^tde, was 
found to be fairly insensitive to the power-law slopes 7 and 
P, but la changes in rb or e can change A^'tde by a factor 
of a few. The sensitivity to rb really reflects a sensitivity 
to the ratio «k/t'osc, as the wandering time can jump by 
~ 1 — 2 orders of magnitude when v^/vcsc rises above a 
value ~ 0.5 — 0.6 (where Ueac here is the escape velocity of 
the stellar bulge). This sensitivity to kick velocity is due to 
the decreased effectiveness of dynamical friction once the 
SMBH begins passing through the galactic center at high 
velocities, giving the black hole more time to disrupt bound 
cloud stars at an observably offset distance. 

Figures [5] and |6] illustrates how our results vary with 
assumptions about the kick velocity distribution, final spin 
amplitudes of the SMBHs, and existence of super-Eddington 
outfiows. In both figures, bound cloud disruptions are repre- 
sented as thick lines and unbound stellar disruptions as thin 
lines. The total number of disruptions is shown as a solid 
line, while those with an observable spatial offset are shown 
with a dotted line, and those with an observable kinematic 
offset are shown with a dashed line. Unless otherwise noted, 
discussion of TDE rates in this section refers only to SMBHs 
which remain bound to their host galaxy. 

Figure [S] displays dP/dluMBH, the number of TDEs ob- 
served by LSST per year per logarithmic black hole mass, 
for our models without super-Eddington emission. Both the 
unaligned (< 180°) and moderately aligned (< 30°) progen- 
itor spin models produce interesting values of P. Both bound 
and unbound disruption rates are dominated by the highest 
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Figure 3. The number of stars, A^tdEi disrupted per year for the 
SMBHs in previous figures. The thick lines refer to disruptions 
from the bound cloud, while the thin hnes refer to disruptions 
of unbound stars. Colors and line types represent the same black 
holes as in Figure [Tl 
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Figure 4. The rate of unbound stellar disruption, A^tdEi as a 
function of radius for the SMBHs in previous figures. 



black hole mass permitted by its spin amplitude to disrupt 
stars; for the first row (fiducial a values), this corresponds 
to while for the second row (a = 0.95) it is 10** 

and for the third row (a — 0.0) it is IO^'Mq. Almost all 
bound cloud disruptions have an observable spatial offset, 
while unbound disruptions never have an observable offset. 
On the other hand, a higher fraction of unbound disrup- 



tions possess an observable kinematic offset relative to the 
bound cloud disruptions. Both these correlations are eas- 
ily explainable: due to high orbital eccentricity, the SMBHs 
in our sample spend the majority of their time far from the 
galactic nucleus, so most bound cloud disruptions occur with 
a large physical offset and low velocity. At the same time, 
virtually all (see Figure |3} unbound disruptions occur dur- 
ing perigalacticon, where the SMBHs move at their highest 
velocities. The highly aligned (< 10°) progenitor spin model 
produces a negligible number of disruptions; high kick veloc- 
ities are suppressed, and SMBHs escape the galactic nucleus 
too infrequently to disrupt significant numbers of stars. 

In Figure |6l we display dF/dluMBH for models with 
super-Eddington fiares. The results are similar to those in 
FigureO although dF/dlnMeH is everywhere greater than or 
equal to its values in the previous figure. Two special points 
of contrast are the large increase in observable disruptions 
at the low end of the SMBH mass function, and the (corre- 
sponding) increase in disruptions for the highly aligned pro- 
genitor spin model. The addition of super-Eddington flares 
has, as expected, little effect on values of dr/dlni\fBH above 
IO^Mq, but disruption flares become dramatically more vis- 
ible for 10^ M0 and SMBHs. 

Figure [7] displays dFs/dlnMBH, the mass dependence 
of the total observable (spatially offset) rate Fg, integrated 
over all kick velocities and all galaxies in our mock cata- 
log, and given flducial spin values. When we integrate over 
black hole mass, we flnd that two of our kick velocity dis- 
tributions produce a robustly observable (~ 10) number 
of disruptions per year assuming super-Eddington flares, 
while the third produces a marginal number of TDEs, of 
order unity per year. Likewise, progenitor spin distributions 
aligned to within 180° or within 30° produce ~ 1 flare per 
year with an observable spatial offset if we are only able to 
observe disk emission. Higher-mass SMBHs contribute the 
most to observable disk fiares, due to the lower temperatures 
and higher optical luminosities of their disks, while super- 
Eddington accretion fiares are dominated by the lower-mass 
end of the SMBH distribution. Although the rate enhance- 
ment from inclusion of super-Eddington outfiows is almost 
a factor of 10, this is considerably lower than the compara- 
ble factor in SQ09. The reason for this disparity is that the 
brightest super-Eddington outfiows correspond to the deep- 
est plunges (lowest Vp) into the tidal disruption region. SQ09 
considered a logarithmically fiat distribution of Vp, while we 
took a constant Tp = rt, for the reasons explained in §3. 

A variety of observable TDE rates are displayed in Ta- 
ble 1. These numbers have been integrated over galaxy type, 
kick velocity distribution, inclination angle, and black hole 
mass function, and indicate that the ultimate observability 
of recoil-induced TDEs will depend strongly on both the ex- 
istence of super-Eddington fiares, and the average distribu- 
tion of pre-merger spin alignments. In this table, F values for 
fiducial SMBH spins are shown, with the a = and a = 0.95 
cases appearing as lower and upper limits in parentheses. It 
is only in the case where super-Eddington fiares do not ex- 
ist and substantial progenitor spin alignment occurs where 
we expect LSST to observe negligible numbers of spatially 
offset TDEs per year. We note that if the progenitor spins 
are unaligned, or even aligned with scatter ^ 30° , the tidal 
disruption rate from recoiling black holes is almost 1% of the 
total TDE rate for all galaxies. For most of our models, the 
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Figure 5. Mass dependence of total TDE rates, F, for a variety of kiclc distributions and mean spins. In the left (green) column, we 
show the unaligned spin case; in the middle (blue), spins aligned to within 30 degrees; and the right (red), spins aligned to within 10 
degrees. The top row has fiducial final spin amplitudes, while the middle has a = 0.95 and the bottom row has a = 0.0. Within each 
frame, the thick lines represent disruptions of bound stars while the thin lines represent disruptions of unbound stars. The solid lines 
represent total number of disruptions, while the dotted lines represent disruptions with an observable spatial offset and the dashed lines 
represent disruptions with an observable kinematic offset. In this plot only disk emission is considered. 



number of kinematically offset TDEs, Fk, is comparable to 
Fs, although we note again that the theoretical basis for ex- 
pecting appropriate absorption lines in a super-Eddington 
outflow is less secure than for a simple spatial offset. We 
have also included in Table 1 the rates of spatially and kine- 
matically offset TDEs for SMBHs which escape their host 
galaxy altogether, labeling these as Fg^eac and Fk.esc- Only 



in the case of unaligned spins and super-Eddington outflows 
are Fs,esc ~ Es and Fs,esc ~ E^; in all other scenarios the 
number of observable TDEs due to ejected SMBIfs is at 
least a factor of 7 smaller than the number due to bound 
SMBHs. 

For most of the models we have considered, a large ma- 
jority of the TDEs associated with recoiling SMBIfs occur 
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Figure 6. The same as the previous figure, but assuming the existence of super-Eddington outflows. 



for black holes bound to their host galaxy. This is due to two 
factors: both the relatively low fraction of SMBHs recoiled at 
escape velocity (see Figure [8] for a plot of the SMBH escape 
fraction, /esc), and the smaller, more rapidly decaying bound 
clouds of those low-mass SMBHs which do escape. This high- 
lights the importance of searching for SMBHs bound to the 
bulge or halo of their host galaxy; although the intergalactic 
TDEs of the KM08 scenario offer a cleaner signal, they are 
intrinsically much fewer in number. 

Finally, it is worth summarizing the primary assump- 
tions we made in this work, where we have tried to err 
on the side of conservatism. To simplify our calculations. 



we neglected emission from the unbound TDE debris, al- 
though that can substantially increase optical (non-super- 
Eddington) emission for low-mass SMBHs (SQ09). Our 
SMBH wandering lifetimes were likely reduced by the fact 
that we limited ourselves to axisymmetric stellar bulge ge- 
ometries, and even more importantly only considered spher- 
ical dark matter haloes. Our simple choice of stellar mass 
function is slightly conservative for calculations of TDE rate. 
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Table 1. Total number of offset TDEs observable by LSST per year (Fs, Fjj) in different models. The fiducial 
spin case is given first, and is bracketed in parentheses by the a = and a = 0.95 cases. 



Spin Alignment [°] 


super-Eddington? 


TslTDEs yr-l] 


TklTDEs yr-i] 


rs,esc[TDEs yr-i] 


rk,csc[TDEs yr-i] 


10 


no 


(0, 0) 


(0, 0) 


(0, 0) 


(0, 0) 


10 


yes 


0.7 (0.7, 0.7) 


0.5 (0.5, 0.5) 


(0, 0) 


(0, 0) 


30 


no 


0.7 (0.3, 1.4) 


1.7 (0.2, 5.6) 


(0, 0) 


(0, 0) 


30 


yes 


9.8 (9.5, 10) 


5.9 (4.4, 9.7) 


1.0 (1.0, 1.0) 


0.8 (0.8, 0.8) 


180 


no 


2.8 (0.4, 6.4) 


3.4 (0.3, 20) 


0.39 (0.3, 0.5) 


0.28 (0.2, 0.3) 


180 


yes 


13 (11, 16) 


7.6 (4.5, 24) 


5.6 (5.6, 5.6) 


4.3 (4.3, 4.3) 
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Figure 7. Mass dependence of the total, galaxy- and velocity- 
averaged rate of spatially offset TDEs, Ts. 180°, 30°, and 10° 
progenitor spin alignment correspond to solid green lines, dotted 
blue lines, and dashed red lines, respectively. Thick lines represent 
disk emission only, while the thin lines correspond to disk plus 
super-Eddington outflows. 

5 CONCLUSIONS 

In this paper, we have demonstrated that super-Eddington 
flares from recoiling black holes, if they exist along the lines 
envisioned in SQ09, will be observably offset to LSST in 
numbers ranging from ~ 1 to ~ 10 TDEs per year. This 
is true for a broad range of assumptions about kick ve- 
locity distributions and galactic structure. This subset of 
transients would contain important scientific value as evi- 
dence of black hole recoil, and could potentially constrain 
the Wk distribution. If super-Eddington flares do not exist 
or if they differ significantly from the SQ09 picture (for ex- 
ample, if /out ^ 0.1), optical emission from the accretion 
disks of TDEs around recoiling black holes will still be ac- 
cessible to LSST, although here the case is more marginal. 
If a large fraction of local universe SMBH mergers proceed 
without significant spin alignment, the prospects for opti- 
cal detection of disk emission from recoiling TDEs are rela- 
tively good, but moderate amounts of alignment would likely 
suppress this. Importantly, the majority of recoiled SMBHs 
will remain bound to their host galaxies, making photomet- 
ric subtraction critical for identification of recoil-associated 
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Figure 8. Fraction of recoiled SMBHs which escape into inter- 
galactic space, as a function of black hole mass, for the three 
different kick velocity distributions. As in the previous figure, the 
green solid line represents the 180 degree alignment distribution, 
the dotted blue the 30 degree, and the dashed red (not visible; a 
negligible fraction of SMBHs from this distribution escaped their 
host galaxies) the 10 degree. 



disruption flares. Depending on the nature of the super- 
Eddington phase of accretion, a comparably large popula- 
tion of kinematically offset flares is potentially detectable, 
but would require spectroscopic foUowup to be realized. 

We have also shown that confusion with TDEs from sta- 
tionary SMBHs will not be a major challenge in the detec- 
tion of off-nuclear TDEs, leaving supernova contamination 
as the main concern. If the scientiflc potential of spatially 
offset TDEs is to be utilized, it will be necessary to con- 
struct transient survey pipelines which do not employ the 
typical "galactic center" cut when searching for TDEs. Al- 
though the challenges inherent to TDE identification have 
been discussed elsewhere (van Velzen et al 2010), the distinc- 
tive lightcurve and color evolution of tidal disruption flares 
are helpful in separating them. The large number of TDEs 
expected to be observed by time domain surveys in the com- 
ing decade will calibrate our understanding of these events, 
so that once LSST is online, it may be able to confirm the 
SMBH recoil predictions of numerical relativity. 
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